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ABSTRACT 

I 

\i' 

■An  algorithm  is  proposed  for  solving  one- dimensional  free  boundary 
problems  with  change  of  phase .  The  technique  consists  of  solving  the  heat 
equation  in  progressively  increasing  rectangles  whose  size  is  controlled  by 
the  Stefan  condition,  convergence  of  the  scheme  is  shown  and  an  estimate  of 
the  rate  of  convergence  is  given . 
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SIGNIFICANCE  AND  EXPLANATION 


Stefan  type  problems  arise  as  descriptions  of  phenomena  such  as  melting 
of  metals,  solidification  of  alloys,  crystal  growth,  permafrost  behavior.  Of 
particular  importance  in  such  problems  is  the  shape  and  evolution  of  the  free 
boundary  or  interface . 

This  paper  suggests  a  method  of  constructing  the  free  boundary  by  solving 
the  heat  equation  in  a  sequence  of  increasing  rectangles .  The  interface  is 
then  approximated  by  piecewise  vertical  segments. 

The  simple  geometry,  and  boundary  conditions  suggested,  can  be  used  to 
perform  efficient  numerical  calculations. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


AN  ALGORITHM  FOR  THE  ONE-PHASE  STEFAN  PROBLEM 
E.  DlBenedetto^  and  R.  Spigler*2^ 

1 .  Introduction ; 

Consider  the  following  one-phase  one- dimensional  Stefan  problem 
Lu  =  uxx  -  ut  “  0  in  Oj  =  {0  <  x  <  s(t)}  *  (0,T] 

u(x,0)  =  h(x) ,  0  <  x  <  a (0 )  -  b,  (b  >  0) 

(SPi  ux(0,t)  -  g(u(0,t) ,t) ,  0  <  t  <  T 

u(s(t),t)  =0,  0  <  t  <  T 

ux(s(t),t)  -  -s(t),  0  <  t  <  T 

where  x  *  h(x),  (C,t)  +  g(€,t)  are  given  functions  on  f 0 , bj  and  a  *  [0,T] 
respectively . 

Under  suitable  assumptions  on  h(*)  and  g,  ( SP )  admits  a  unique  classical 
solution.  For  such  results  we  refer  to  the  survey  article  [13]  and  other  papers  given 
in  the  extensive  bibliography. 

The  aim  of  this  paper  is  to  propose  an  algorithm  to  construct  the  solution,  which 
consists  in  solving  the  heat  equation  in  progressively  increasing  rectangles,  whose 
si*e  is  controlled  by  the  Stefan  condition  ux(s(t),t)  “  -s(t). 

Such  an  algorithm  arises  as  a  natural  modification  of  Huber's  method  [15,10,1] 
and  can  be  described  in  a  simple  fashion  as  follows. 

First  the  interval  [0,T]  is  divided  in  n  intervals  of  length  0  *  T/n,  then 
for  t  e  [0,9]  we  set  Sg(t)  “  b  and  solve  the  problem 
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uxx  "  ut  ■*  0  in  =  {O  <  x  <  b)  »  (0,6j 
u’(0,t)  -  g(u1(0,t),t),  0  <  t  <  8 

(P.,)  u'fX/O)  «  h(x),  0  <  x  <  b 

u^b.t)  -  0,  0  <  t  <  8  . 

Ms  compute  the  number  u\sg (9 )  -  b,8)  and  determine  the  rectangle 
Sj  Mo  <  X  <  x2  -  b  -  u’(b,8)6}  X  (9,28]  , 

setting  Sg(t)  “  x^  for  t  e  (9,28].  in  Rj  we  solve  a  problem  similar  to  (P,),  and 
proceed  in  this  fashion . 

The  convergence  of  schemes  where  at  each  time  step  the  free  boundary  is 
approximated  by  a  vertical  segment  was  conjectured  by  Datzeff  [5,6]  . 

A  proof  of  convergence  has  been  given  by  Fasano-Primicerio-Fontanella  [11]. 

Their  scheme  however  is  somewhat  more  complicated  than  the  one  we  propose  here,  both 
in  the  construction  of  the  sequence  of  rectangles  and  in  the  boundary  conditions  on 
xj>(j  "  1)9  <  t  <  j8  which  are  not  homogeneous,  being  given  as  a  relationship  linking 
the  distance  t  -  (j  -  1)8  with  the  values  ux-1 <xj_i » ( j  "  D®>* 

Thus  the  scheme  we  hove  described  has  a  two-fold  simplicity t  the  rectangular 
geometry  and  the  homogeneity  of  the  boundary  data  on  the  approximating  free  boundary. 

<te  treat  the  problem  for  boundary  data  on  x  -  0  of  variational  type  since  such 
a  condition  is  the  'natural"  one,  as  pointed  out  and  discussed  in  [7] .  The  method 
could  handle  the  nlrichlet  boundary  data  as  well. 

He  give  an  estimate  of  the  speed  of  convergence  which  turns  out  to  be  of  the 
order  of 

As  a  related  work  we  mention  the  methods  of  [14]  based  on  enthalpy 
considerations,  which  yield  a  slightly  poorer  rate  of  convergence,  of  the  order  of 
[8  in  8”1]  1y2  . 

The  methods  of  proof  are  simple  in  that  we  exploit  both  the  rectangular  geometry 
and  the  homogeneity  of  the  data  to  represent  the  approximating  solutions  by  means  of 
elementary  heat  potentials. 
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Section  2  contains  the  precise  description  of  the  algorithm,  assumptions  and 
statement  of  results.  In  Sections  3  and  4  we  produce  basic  estimates  and  prove  the 
convergence  of  the  approximating  solution,  to  the  solution  of  (SP) . 

The  error  estimate  is  given  in  Section  S,  He  conclude  the  paper  by  discussing 
some  variants  of  the  scheme. 


2.  Assumptions  and  statsasnt  of  results! 


Throughout  tha  paper  we  will  make  tha  following  assumptions  on  tha  data. 

[A,]  x  ♦  h(x)  ia  a  poaitiva  Lipschitz  continuous  function  on  [0,b],  with 
Lipschitz  constant  H,  and  h(b)  «  0. 

Ia2]  (£,t)  ♦  g(€,t)  is  non-positive  on  S  x  (0,T],  continuous  with  raspact  to 
t  e  10, T) ,  Lipschitz  continuous  in  £,  uniformly  in  t,  with  Lipschitz 
constant  Gj  and  g(h(0),0)  •  h’(0).  Moreover  thara  exists  a  non-negative 
constant  G2  such  that 

lg(C,t)  I  <  Gn  151  +  a2,  (C,t)  e  *  *  (0,t1  . 

for  n  *  1,2,...,  sat  9  «  T/n  and  cons ider  tha  sequence  of  problems  Py 


1 


1,2, ...,n,  defined  by 
Lu^  »  0  in 

u3(x,(j  -  1)6) 


Rj  =  {0  <  x  <  x^)  x  {{j  -  1)0  <  t  <  jQ} 
fu^_1(x,(j  -  1)6),  0  <  x  <  x^_1 


'V 


;o  • 

u£(0,t)  •  g(u^ (0,t) ,t) , 

u^(x^,t)  -  0, 


*j-1  <  X  <  X) 

(J  -  1)6  <  t  <  j 0 
(J  -  1)0  <  t  <  j6 


where  tha  sequence  (x^)™^  is  recursively  defined  by 

Xq  -  x,  »  bf  Xj  -  x^_,  -  u^”1(Xj.1,(J  -  1)0)*9,  i  -  2,3,...,n 

and 


u°(x,0)  -  u’lx.O)  5  h(x)  in  [0,b]  . 

By  virtue  of  [A^l  -  [A2I  each  (Pj)  admits  recursively  a  unique  classical 
solution  u^  whose  derivative  uj  exists  up  to  the  lateral  boundaries  of  R j . 
Conn«.juently  the  sequence  (x^)  Is  well  defined. 

Setting 

Sg(0)«  b»  Sg(t)  «  Xj  for  (j  -  1)0  <  t  <  j0,  j  «  1,2, ...,n  , 
we  obtain  a  right-continuous,  piecewise  oonstant  function  defined  in  [0,T]  .  By 
elementary  considerations  and  the  maximum  principle  [7,8]  u^(x,t)  >  0  and  the 
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3.  Some  basic  estimates! 

Let 

r<x,t>5,T)  »  — - - - - exp 

2/*(t  -  T) 

be  the  fundamental  solution  of  the  heat  equation  and  let  G(x,tj?,T),  N(x,t}£,i)  be 
the  Green’s  and  Neumann's  functions  respectively,  defined  by 


G(x,t»C,T)  »  r<x,tj5,T)  -  r(-x,t>C,T)  , 

N{x,t>5,f)  -  r<X/tie,T)  +  rt-x,tiC,T)  . 

In  the  rectangle  R  ^  the  solution  u3  of  (P^ )  can  he  implicitly  represented  as 


(3.1)  uj(x,t)  =  /  N(x,t,C,(j  -  1)6)uj(C,(j  -  1)8)dC  - 

0 


V  , 

-  /  N(X,tj0,t)g(u3(0#X) ,T)dT  + 

(j-1)8 


, 

♦  /  N(x,tix  ,t)u'(x.,t)dt 
( j-1 )8  3  x  3 


for  (x,t)  e  Rj .  Taking  the  derivative  with  respect  to  x  in  (3.1)  and  letting 
x  ♦  Xj  we  obtain 

Xj 

(3.2)  ~  u^  (x . ,t)  »  /  G(x.,tiC,(j  -  1)8)u^(?,(j  -  1 )  8 ) )  dt  - 

*  *  J  Q  J  X 

t  , 

-  /  N  (X  ,t|0,T)g(u3((),T)  ,T)dT  + 

(j-1 )9  x  3 
t  . 

+  /  N  (X,,t>X,,T)U3(x,,T)dT,  (j  -  1)8  <  t  <  j8  . 

(j-1 >6  x  3  1  x  3 

The  calculations  leading  to  (3.1)  —  (3.2)  are  routine  and  we  refer  to  [4,9,12]  for 


details . 


f 


Let  us  fix  1  <  j  <  n  and  (x,t)  *  Rj  and  integrate  the  Green  Identity 

3  3 

(Nu?  -  uN?)  -  npj  (Nu)  *  0 


over  R.  ,  1  <  k  <  j.  Since  in  R. 


V  we  are  away  from  the  singularity  we  obtain 


H 

H 


xk  ** 

(3.3)  /  N<x,tte,k9)uk(5,k6)d5  -  /  N (x ,t 1 5 , (k  -  1)9)uk(C,(k  -  l)9)d£ 

0  0 


>t0  k9 

-  -  J  N(x,tj0,T)g(ua(0,T),T)dT  +  /  N(X,t>X,./T)uk(X  ,T)dT 

(k-1 )9  8  (k-1)9  k  *  k 


for  k  -  1,2,  -  1) . 

By  virtue  of  our  definition  of  (Pj),  the  second  integral  on  the  left  hand  side  >f 
(3.3)  can  be  rewritten  as 


I 

\ 


Vi 

-  J  N(x,ti£,(k  -  1  }9)uk_1  (£,  (k  -  1)9)d£  . 
0 


Consequently  adding  the  identities  (3.3)  for  k  -  1,2,...,(j  -  1)  with  (3.1) 
recalling  the  definition  of  ue(x,t)  we  obtain 


and 


b  t 

(3.4)  U0(x,t)  -  /  N(x,tjC,0)h(S)d£  -  /  N(x,tlO,T)g(U.(0,T),T)dT 

0  0  8 


j-1  k9 

*11  N(x,t»x  ,T)uk(x  ,T)dT  + 

k-1  (k-1 ) 9  *  *  * 


♦  /  N(x,t)X  ,T)u^(x.,T)dT 

(1-1)9  3  x  J 


< j-1 )9 

Lemma  3.1;  For  each  9  the  following  estimates  are  valid 


(a) 


2G  ft  2G./T//i 

'u9(0,*)It  <  [Hb  +  -f—  ]e  SC  . 

/w  u 


(b> 


*u9 '  «  Co  * 

m,s 
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f 


i 

P 


Proof  i  By  virtue  of  the  maximum  principle  u-  >  0  in  V  and  uj(x4,t)  <  0, 

B  8g  3 

t  e  ( < J  -  1)8,  j8],  therefore  dropping  the  non-positive  terms  on  the  right  hand  side 
of  (3.4)  and  letting  x  *  0  we  obtain 


b  t  g(u0(O,T),T) 

<  ug(0,t)  <  /  2r(0,t>5,0)h(C)d£  -  —  /  - ■-  - -  dt 

0  /»  0  /t  -  T 

Hh  ,b  1  3  G1  ,fc  ua(°*T>  _ 

<  —  /  —  exp[-C2/4t] d$  +  —  /  ■  ■  —  dT  + -  /t  . 

/»  o  /t  /a  o  /t  -  t  /» 


Statement  (a)  is  now  a  consequence  of  standard  calculations  and  Gronwall's 
inequality. 

Statement  (b)  follows  from  the  maximum  principle  applied  recursively  to  (P^). 
Lemma  3.2;  For  each  j  -  1,2,...,n 

(a)  |u£(x^,t)|  <  HeT  +  BGfe*  SC,,  f  -  T//*b2 

3 

(b)  |-j^ue(x,t)|  <0,,  (x,t)  e  Vg  , 

..  8 

where  G  -  G,C0  +  Gj,  and  H  -  max{H,G>  . 

Proof:  We  employ  an  induction  argument  by  making  use  of  formulae  (3.2).  First  we 
prove  that  if  for  some  j  -  1,2,...,(n  -  1)  we  have 


l“^(x,30) I  <  P,  0  <  x  <  Xj 


for  some  positive  constant  P,  then 


(3.5)  |uj> 


^+1(x1+1,t)|  <  Pe9//?b  +  8G  -2—  e6/^b  ,  t  e  ( j8,( j  +  1) 


max  lg(ufl(0,t) ,t> |  <  G.C.  *  G,  :  G 


Consider  (3.2)  written  for  the  Integer  j  +  1 

0.6)  *3 

and  estimate  the  integrals  1^,  i  ■  1,2,3  separately  as  follows. 


V.  — sgi 


1+1 


1,1  »  I  /  G(x^,tj£,j0)u3  1(C,j6)dC|  -  by  definition  of  (P^ ) 

0  3 

Xj  xj 

-  | /  G(x  ,ti?,j0)u^(£,j0)d£|  <P  /  |G(x  ,t;S,jS)|d€ 
0  3  0  3 

*  2 

.  p  f  -n  p 
<  —  J  e  dn  -  -  . 

/»  0 


To  estimate  Ij»  1-j  we  recall  the  following  elementary  estimates  on  Nx 


lVx  t,0,T)|  < 


/Jb: 


2  ' 


(xj+1  >  b  >  0) 


*  3  +  1  '  J  +  1  2/^b2 


Therefore 


Putting  together  these  estimates  as  parts  of  (3.6)  we  have 


for  all  t  e  < 3 9  ,  (  j  ♦  1)6). 

Consequently  by  Gronwall's  inequality  (3.5)  follows  at  once. 

Consider  now  the  problem  (Pj)  j  =  1.  since  x  ♦  h(x)  is  Lipschitz  continuous 

in  (0,b)  ,  h'(x)  exists  for  a  .e  .  x  e  [ 0 , b]  and 

ess  sup  |h'(x)|  <  H  <  H  »  max  {h,g}  . 

1 0 ,  b] 

Therefore  by  the  previous  argument 
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Iui*x1't)  I  *  "e°  +  ®a«a»  t  e  (0,8) 

where  for  simplicity  of  notation  we  have  set  a  «  0//»b2.  The  function 
(x,t)  ♦  u^(x,t)  s  v(x,t),  will  satisfy  the  Dlrichlet  problem 
Lv  »  0  in  R, 

v(0,t)  «  g(u0(O,t),t)>  I  g(u0(O,t)  ,t)  |  <  G,  te  (0,8) 

|v(x.,,t)|  <  He°  +  8G0ea,  t  8  (0,8) 

v(x,0)  -  h'(x),  |h'(x)l  <  H  s.e.  x  8  (0,b)  . 

Consequently  by  the  maximum  principle 

|uNx,8)|  <  He°  +  BGoea 


and  by  (3.5) 


|u^(x2,t)|  <  He20  *  8Gae2a  8Gae° 

<  He2®  +  0G(2«)e28,  t(6,28)  . 


Proceeding  In  this  fashion  we  obtain 

|u£(xj,t)|  <  He^°  +  8G(j«)e^a  . 

Now  ja  »  -  -  j8  <  T,  and  therefore  the  lemma  is  proved.  Next  we  Introduce  the 

/xb2 

function  t  *  a0(t)  defined  by 

(3.7)  S  (t)  -  x.  -  u^(x  ,J8)  (t  -  (j  -  1)8),  te  II)  -  1)8, j8)  . 

0  J  x  j 

For  t  »  j8,  s0(t)  ■  xj+i'  8g(0)  b»  80  that  the  of  a0(*)  iB  obtained  by 

connecting  the  points  (b,0 ) , (x. ,8) , . .. , (x  ,{n  -  1)8)  for  t  8  (0,(n  -  1)8),  and  by 

2  n 

the  graph  of  (3.7)  for  te  [  (n  -  1)8,t)  .  The  points  (Xj,(j  -  1)8),  j  -  1,2,...,n, 
are  the  lower  vertices  at  the  right  side  of  the  R^'s. 

By  Lemma  3.2 

b  <  Sg(t)  <  b  +  CjT,  t  e  I0,T) 

therefore  the  sequence  (sgt*))  is  egulbounded  and  equilipschitz,  so  that  by  Asaoll- 
Arzela  theorem  a  subsequence  relabeled  with  8,  converges  uniformly  to  some  non- 
decreaainq,  Lipschitz  continuous  curve  t  *  s*(t),  with  Llpschitz  constant  bounded 
by  C,  .  Since  ls0  -  s^l  <  ci®’  *1*°  s0(t)  converges  uniformly  to  s*(t)  . 
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Let  Dg«  be  the  domain  defined  by 

Dg#  5  (0  <  x  <  s*(t)}  *  (0,T]  , 
and  let  u*  be  the  unique  solution  of 

Lu*  “0  in  Dg,  , 

ui(0ft)  -  g(u*(0,t) ,t),  t  e  (0,Tl  , 

(P) 

u*(xf0)  -  h(x),  x  e  (0,b]  , 

u*(s(t),t)  “0,  t  e  [0,t]  . 

We  will  show  that  ufl(x,t)  ♦  u*(x,t)  uniformly  in  S  and  that  the  pair  (u*,s*)  so 
obtained  is  actually  the  unique  solution  of  <SP)  in  the  introduction. 

We  remark  that  as  a  consequence,  in  view  of  the  uniqueness  for  (SP) ,  the 
selection  of  subsequences  is  superfluous. 

The  following  lemma  will  be  needed. 

Lemma  3.3;  Let  (x,t)  *  v(x,t)  be  the  unique  solution  of 

Lv  »  0  in  Og*  > 

vx«>.t)  «  g(u0(O,t)  ,t) ,  te  (0,T)  , 

v(x,0)  »  h(x),  x  e  [0,b]  , 

v(s* <t) ,t)  -  o,  t  e  (o ,t]  . 

Then 

0  <  v(x,t)  <  C^sMt)  -  x) ,  (x,t)  e  Dg*  . 

Proof s  The  lemma  is  proved  by  standard  barrier  techniques  and  the  maximum  principle 
MOI  • 
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4.  Convergence  of  the  sch— a  i 

Lemma  4.1 1  u0(x,t)  ♦  u*(x,t)  uniformly  in  S,  as  8*0. 

Proof i  By  the  triangle  inequality 

lua  *  u*i  4  LI  +  lw  I 
8  -,8  1  »,S  2  «,3 

where  w,(x,t)  “  u0(x,t)  -  v(x,t)  and  w2(x,t)  ■  v(x,t)  -  u*(x,t),  and  v  is 
defined  in  Leana  3.3. 
set 

o(t)  -  min{sg(t) ,s#(t) )i  Bit)  ”  aax{sg(t),s*(t>}  , 

«(t)  -  B(t)  -  a(t>,  t  e  [0,T]  . 

Me  already  know  that  S(t)  *  0  uniformly  in  [0,T] ,  as  6*0. 

Consider  the  rectangle  *j,  1  <  )  <  n.  Me  claim  that  if 
!*,<*,(!  “  1)8)1  <  4^181^  1)8,  x  e  <0,->,  then 

|w,(x,t)|  4  c  l«l  ,  <x,t)  8  (0,«)  X  < ( j  -  1)8, j8)  . 

If  for  t  e  t(j  -  1)9,19],  x^  <  s*(t)  then  w,  solves  the  problem 
Lw,  -  0  in  {0  <  x  <  x^ }  x  ((j  -  1)8, j8]  , 
w,  ( 0 ,t)  -  0,  t  e  ((j  -  1)8, j8]  , 

x 

!w1(x,(j  -  1)8)|  <  c1,5,(j_1)0*  xe  * 

w,(x^,t)  -  -v(Xj,t),  te  ((j  -  1)8, j 0 J  , 

Hence  |w,(x,t)|  <  max(c  max  v(x  ,t)J,  (x,t)  e  Kv  By  Leimna  3 .3  we 

<3  1(1-1)8,381  3  3 

obtain 

lw,(x,t)|  4  (x,t)  e  10,-)  *  t(j  -  1)8, j6]  . 

If  for  t  e  [(j  -  1)8,18],  x^  >  s*(t) ,  then  w,  solves  the  problem 

Lw,  -  0  in  (0  <  x  <  s*(t)}  x  ((j  -  1)6,10]  , 

w,x(0,t)  -  0,  te  II)  -  1)8,10]  , 

l«,(x,(l  -  1)8)1  4  C,l5l(j_1)9,  0  <  x  <  s*[(l  -  1)8]  , 

w,(s*(t),t)  -  u0(s*(t),t),  t  e  ((1  -  1)8,18]  . 

By  Lemma  3 .2  we  have  0  4  u0(x,t)  4  C,(x^  -  x),  (x,t>  e  R^  and  therefore 
|w,(x,t)|  4  C}l«lt,  <x,t)  e  (0,«)  x  ((j  -  1)8,18]  . 
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If  t*  e  ((j  -  1)6, J0)  such  that  x^  »  s*{t*),  than  we  repeat  analogous  arguswnts 
in  the  domains  so  determined. 

Now  since  for  t  -  0,  w1(x,0)  -  0,  x  e  (0,*),  an  inductive  argwent  gives 
|w^ (x,t)  |  <  c1i«it,  xe  (0,-),  te  to ,t]  . 

As  for  Wj,  since  it  solves  the  problem 

Im2  -  0  in  Ds,  , 

w2x(0,t)  ■  g(ug(0,t),t)  -  g(u*(0,t),t)  , 
w2(x,0)  -  0,  x  e  (0,b) 

w2(s*(t),t)  -  o,  t  e  (o ,t] 

it  can  be  dominated  [2,3]  by  the  function 


w  (x,t)  -  /  N(x,t;0,T)|g(u8(0,T),r)  -  g(u*(0,T),T)|dT 

0 


the  unique  solution  of 


U»2  «  0  in  S  , 

"2x(0,t}  ”  #t)  -  g(u*(0,t),t)|,  te  (0 ,T]  , 


w2(x,o)  -  o,  x  e  (o,-)  . 


|w2(x,t)|  <  G1  /  N(x,t»0,T)|u9(0,T)  -u*(0,T)|dT 

0 


deduce  that 


t  |Ug(0,T)  -  U* ( 0 , T ) | 

|uA(x,t)  -  u*(x,t)  |  <  |w.(x,t)  I  +  G.  —  /  - - -  dt 

8  1  1  fi  0  /t  -  T 


sup  |Ug(x,T)  -  U ( X , T ) | 

<  C  .5.  +!l/  - 

^  ^  /*  0  /t  “  T 


i 


And  by  Gr on wall 1 •  Inequality 


(4.1) 


■up  |u8(x,t)  -  U(Xjt) |  <  C  141  e 
x«[0,«) 


-1  G 

✓? 


C2,4,t 


for  all  t  e  [0,T]  .  The  lemma  is  proved. 

Le— a  4.2)  The  pair  (u* ,s* )  coincides  with  the  unique  solution  of  (SP)  • 

Proof i  The  only  thing  that  remains  to  be  proved  is  the  Stefan  condition 

u(s(t),t)  -  -s(t),  t  e  (0,T)  .  Such  a  condition  has  been  shown  to  be  equivalent  to  the 

Integral  identity  [2,  page  85] 


(4.2) 


t  b  s(t> 

s(t)  «  b  -  J  g(u(0 ,t) ,t) dt  ♦  /  h(x)dx  -  /  u(x,t)dx 
0  0  0 


and  hence  it  will  be  sufficient  to  prove  that  a*,  u*,  satisfy  (4.2). 
Integrating  the  equation  Lu^  -  0  over  Rj  we  obtain 

(J6  3  rJ0  3 

-  J  UJ(x  ,t)dx  m  .  J  g(uJ(0,T),T)dt  + 

(3-1)8  3  (j-1)8 


(4.3) 


3  j 

+  J  uJ(x,(3  -  1)8)dx  -  J  uJ(xfj8)dx  . 
0  0 


Also  for  (x,t)  6  Rp,  1  <  p  <  n,  integrate  Lup  ■  0  over  the  rectangle 

{0<x<x)*((p-  1)8»t)  .  It  gives 
P 

.  x 

t  t  p 

(4.4)  -/  uP(x  ,t)dl  -  -  /  g(uP(0,T),T)dT  +  /  uP(x, (p  -  1 ) 9) dx 

(p-1)9  X  P  <p— 1 ) 0  0 


-  /  up(x,t)dx  . 
0 


By  the  definition  of  (P^)  we  have 


:^K 


I 


L  ■' 


M 

4 


therefore  adding  the  identities  (4.3)  for  j  -  1,2,...,p  -  1  and  (4.4)  we  obtain 


(4.5) 


p-1  j0  t 

I  J  U3(X  ,T)dT  -  J  UP<X  ,T)dT 

l-i  (j-tja  x  3  (p-i )a  x  p 


t  b  p 

■  /  9<u8(0,t)  ,T)dT  ♦  /  h(x)dx  -  /  up(x,t)dx  . 
0  0  0 


Me  rewrite  the  left  hand  aide  of  (4.5)  as  follows! 


P-1  18 


P^1  .J0 


-If  U^(X.,T)dT  -  /  >jP(x  ,T)dT  -  -  If  U^(X  ,j6)dT 

j-1  (1-1)9  x  3  (p-1 )8  *  p  j-1  (J-1 )6  x  3 


t  p-1  1®  .  t 

-  /  up(x  ,p®)dl  -If  ei(T)dt  -  J  ep(T)dt  , 

(p-1  )8  x  p  j-1  (j-1)®  (p-D® 


where 


ei(t)  -  u1 (x. ,t)  -  u1(x4<ifl),  (i  -  1)0  <  t  <  i®  ,  i  -  1,2, ...,n  . 

w  X  i  XX 

Me  observe  that  the  numbers  -u^(x. , j 8)  are  the  slopes  of  the  Lipsc-hitz 

^  J 

continuous  polygonal  t  ♦  Sg(t)  for  (j  -  1)®  <  t  <  j®,  consequently 


p-1  1®  .  t  t  . 

-If  u,<x.,l®)dT  -  /  UP(x  ,p®)dr  -  /  —  8g(T)dT  -  8g(t )  -  b  . 


1-1  (1-1)0  3 


(p-1 )® 


*  p 


Carrying  this  in  (4.5)  gives 


(4.6) 


t  b  p 

Bg(t)  -  b  -  /  g(u0(O,T),T)dT  +  f  h(x)dx  -  /  Ug(x,t)dx  + 


P“1  1®  .  t 

*  I  !  e|(T)dx  +  /  ep(T)dT  . 

1-1  (1-1)0  (P-1>® 


By  virtue  of  lemma  4.1  and  the  uniform  convergence  Sg(t)  ♦  s*(t),  letting 
9*0  in  (4.6)  gives 
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(4.7) 


t  b  .*(t) 

e*(t)  «  b  -  J  g(u*(0,T)  ,T)dT  +  /  h(x)dx  -  /  u*(x,t)dx  + 

0  0  0 


‘'1 


-I 


p-i  je  t 

+  11m  {  I  /  e2(T)dt  +  /  e|(T)dT}  . 

0+0  j-t  (3-1)8  (p-1 )0 


Therefore  the  lemma  will  be  proved  If  we  show  that  the  limit  In  (4.7)  le  zero . 

tn  order  to  estimate  the  «g(M  we  will  need  a  representation  for  u£(Xj,t), 
t  e  <(J  -  1)9, J9] . 

Consider  Identity  (3.3).  By  taking  the  derivative  with  respect  to  x  and 
Integrating  by  parts  the  first  two  Integrals  (the  Identity  Nx  -  -G^  Is  used)  we 
obtain 

*k  .  *k 

(4.8)  /  G(x,t»5,k0)u  (C,k0)dC  -  /  G(x,t»C,(k  -  1)0)u  (E,(k  -  1)0)d? 

0  x  0  * 

kO  kO 

m  ~  /  H  (x,t»0,T)g(ufl(0,t),T)dT  +/  N  (x,t»X^T)u*(*  ,T)dT  . 

(k-1 )6  (k-1 )6  x  k  x  k 

Now  by  the  construction  of  uk,  the  second  integral  In  (4.8)  can  be  rewritten  as 

Vi 

-  /  G(x,t»e,(k  -  1)8)u*  (t,(k  -  1  )8)dC  , 

0  x 


i'4 


therefore  adding  the  identities  in  (4.8)  for  k  -  1,2,...,(j  -  1),  changing  the  sign 
and  evaluating  the  sum  at  x  -  Xj ,  gives 

b  xj-i 

(4.9)  J  G(x  ,t jC.Olh1  (£) d€  —  J  G(x  ,t»5,(J  -  1)0)uJ-1(5,(j  -  1)0)dC 

ft  J  ft  J  X 


-  J 

0 


(j-1)0 


N  (x.,t,0,T)g(Ufl(0,T),T)dT  - 
x  j  e 


j-1  k8 

-l  !  N  ( x . , t; x  »  f  )u*!(x.  ,T) 

k-1  (k-1 )0  3  *  x  x 


dT  . 
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Consider  now  the  representation  (3.2).  By  the  way  u3  has  been  constructed,  the 
first  Integral  on  the  right  hand  side  of  (3.2)  reads 


/  G(x  ,t»?,(j  -  1)6)u*-1(€,(j  -  1)0)d?  . 
0  3 


Finally  adding  (3.2)  and  (4.9)  gives 


,  1  b  t 

(4.10)  -  u3(x  ,t)  -  /  G(x  ,t»5,0)h’(C)dC  -  /  N  (x.,t»0,T)g(u.(0,T),T)(jT  + 

303  o  X  J  8 

j-1  le© 

+  I  /  N„(*.ft|X  ,t)u\*  ,t)dT  + 

k-1  Oe-1  )©  x  3  K  x  k 
t 

*  !  N  <x  ,t>x  ,t)u3(x  ,x)dt,  te  ((j  -  1)6, j6] 

< j-1 )©  3  3  x  3 

The  error  terms  e2( t)  J  “  2,3,...,<p  -  1),  can  now  be  written  as 


^e3(t)  «-j  (u^.t)  -  u^Xj.jflM  -  /  (Gfx^tiC.O)  -  GfXj.jflK.OJIh'fOdC  - 
t 

"  /  (N  (X  ,t»0,T)  -  H  (X  .,je,0,T)]g(u  (0,T),T)dT  + 
ft  *  J  X  J  0 


+  /  N  (x  ,j9,0,T)g(u  (0,T),T)dT 
t  x  3  8 


j-1  k0 

+  I  /  IN  (X  ,t»x  ,t)  -  N  (X  ,j0,x  ,T))uNx.  , 

k-1  (k-1 )0  x  3  k  x  3  k  x  k 


T)dt  + 


/  tN  (x  ,t|X  ,T)  -  N  (X  ,j0|X  ,t))u3(x  ,T)dT 

(j-1)0  *  3  3  XJ  j  Xj 


Jw 

/  Nx(*j,j0|Xj,T)u3(x3,T)dT  - 


J1  +  J2  +  J2  *  J3  +  J4  +  J4  • 


By  the  nean  value  theorem  and  standard  estimates  on  the  exponential  function  we  have 


|G ( xj ,tiC/0)  -  G(xj> J8>C,0) |  <  (K1  +-37j>8  <  Kj8  ♦  X2<i  -  1)'3/28"1/2  . 


Here  and  In  what  follows  Kg  denote  constants  Independent  of  j .  Consequently 

U,!  <  HbO^S  ♦  K2(j  -  1)“3/28“1/2)  . 

As  for  Jj  analogous  arguments  yield  the  estimate 

|N  (x  ,t>0,T)  -  N  (x  ,56,0,1)1  <  K  8  , 

X  J  X  J  £ 

where  we  have  used  the  fact  that  x^  »  b  >  0.  Therefore  J2  can  be  estimated  by 

|J2I  <  GKjT8  . 

Let  us  now  turn  to  J3 .  First  ws  note  that  for  k  ■  -  2), 

In  (x  ,tix,t)  -  n  (x  ,i8»x  ,t)|  <  (j8  -  t)  sup  |n  <x  ,P;x  ,t)|  < 
x  j  x  x  3  x  xt  J  K 


<  K48  k58  sup 


(xj  ~  V 

1  '  *fc  "  4  ( P  -  T) 

r  /a  ®  • 


t<p<j8  (p  -  T) 


5/2 


Now  since  (x^.ke),  (x^,<l  -  1)8)  belong  to  the  maximal  extension  of  the  graph 
<t,8g<t)),  we  have  <x^  -  x^)  <  ZC^IP  -  T),  for  k  -  1,2,...,(j  -  2),  and  hence 

j-2  k8 

U,l  <  l  I  In  <x.,t;x  ,r)  -  N  <x  ,10>x  ,r)||u  <m  ,T)|dr  + 

3  k-1  (k-1 )8  x  3  k  x  3  K  x  * 


<1-1)8 


♦/  |N  (x  ,t)X  fT)  -  N  <x  ,}8>x  ,T)||u3”  <*.,_1»T>|dT 

(j_2)0  *  j  j  1  J  *  J1 


<1-2)8  ,  <1-2)8 

<  C|K  0  J  dT  +  2C,K.0  /  [<}  -  1)8  -  t)  dT  ♦ 

0  0 


<1-1)0 

♦  C,  J  |N  <X  ,t»X,  ,t)  -  N  <X  ,}8,x  ,T)|dT  . 

1  ( j -2 ) 8  x  3  3-1  x  3  3-1 


In  order  to  estimate  the  last  Integral  recall  that 


*x<Xj 'P,Xj-1 »T) I  *  K6<p  ■  T)  '  P  >  T  • 


Substituting  this  in  the  estimate  of  Jj  we  obtain  after  some  algebra 

|J3I  <  K?/3  . 

Analogous  considerations  yield  |J£|  +  I J^l  +  |J^|  <  Kg/0.  Therefore  the  error  terms 
Sg(t)  can  be  estimated  by 

I  <  Kg/e  +  K10(J  -  i)"3/2e",/2  . 

For  the  error  Eg  in  (4.6)  we  have 

p-1  J9  t 

I  Eg |  <  l  I  |el(T>|dT  +  /  |e|(T)|dT  < 

j-1  (j-1)0  (p-1 >9 


9  Hi®  a 

/  |eI(T)|dT  +  I  /  leimldT  . 

0  8  j-2  (j-1  )9 


Wow  trivially  |eft(T)|  <  2Cf .  Therefore 


(Eg |  <  2C  9  +  K  ^  !  8  +  K  /5  f  (j  -  1)‘3/2  . 

8  1  9  j-2  10  j-2 


Since  f  0  <  p6  <  T,  and  the  series  I  (j-1)  3^2  converges, 
j-2  j>2 


obtain 


(4.11) 

This  proves  the  lemma . 


|Eel  <  k,,/5  . 


•  41 


5.  The  error  estimate; 

The  purpose  of  this  section  is  to  give  sn  estimate  of  the  speed  of  convergence  of 
the  approximate  interface  Sg(t)  to  the  true  interface  s(t) .  In  view  of  Lemma  4.1 
and  (4.1)  this  will  also  complete  the  proof  of  the  theorem. 

Lemma  S , 1 r  There  exists  a  constant  C  depending  upon  H,  b,  G,,  Gj,  T  such  that 

,a8  "  *  c/S  * 

Proof;  Subtract  (4.6)  from  (4.2)  and  use  (4.11)  to  obtain 


|s(t)  -  se(t)|  <  /  |g(u(0,T) ,x)  -  g(ue(0,T),t)|dT  + 

0 

s(t )  88{t> 

+  |/  u(x,t)dx  -  /  Ug(x,t)dx|  + 


s6(t) 


+  |  /  Ug(x,t) dx  -  /  Ug(x,t) dx|  +  K  /5  » 

0  0 


where  u0  is  the  unique  solution  of  the  problem 

Lun  «  0  in  V-  3  (O  <  x  <  s.(t)}  *  (0,T]  , 

8  Sg  8 

Ugx(0,t)  -  g(Ug(0,t) ,t) ,  0  <  t  <  T  , 

Ug(x,0)  -  h(x),  x  e  (0,b)  , 

Ug  (s g  (t) ,  t )  -  0,  0  <  t  <  T  . 

On  the  basis  of  the  maximum  principle  and  standard  barriers  estimates  [10,7],  as  in 

Lemma  3.3,  we  have 

0  <  u(x,t)  <  (s(t)  -  x),  (x , t )  e  dt  , 

0  <  u0 ( x , t )  «  C1(Sg(t)  -  x),  (x,t)  e  Pj  . 

8 


Proposition  5.1;  There  exist  constants  B^  and  Bj  depending  upon  H,  b,  G ^ ,  G2 
T ,  such  that 


-  <Jt». 


3 


s(t) 


«8<t) 


|/  u(x,t)dx  -  /  u9(x,t)dx|  <  B  / 

0  0  0 


SUP  |s(P)  -  Sfl(P)  | 


0<P<T 


dr  + 


ft  -  T 


+  B%; 


t  l«l 


dT  . 


Proof  of  proposition  5.1 «  Set 

a(t)  -  min{s(t) ,s9(t)}»  B(t)  -  m«x(s(t) ,s0(t) } , 

5(t)  -  B<t)  -  3(t),  t  e  [0,T]  . 

Obviously  a(»),  ?(•)  are  non-decreasing  Lipschitz  continuous  functions  with 


As  for  Vj,  it  can  ba  represented  explicitly  (9J,  by 


v,(x,t)  -  /  N(x,ti0,T)|g(u(0,T),T)  -  g(u.(0,T)  ,x)  |dT 

0 


Therefore  toy  virtue  of  the  Lipschitz  continuity  of  g(*,t)  and  (4.1)  we  obtain 


Gi  ,4,t 

V2(x,t)  <  —  C,  /  - - dt 

/»  *  0  / 1  -  t 


By  an  argument  of  [*.  page  87/  v?  can  be  dominated  by  z(x,t)  +  z(-x,t)  where  z 
so Ives 

Lz  -  0,  {-a(t>  <  x  <  ■}  x  {O  <  t  <  T>  , 

z(-a(t),t)  -  C^Jft)  ,  0  <  t  <  T  , 

z(x,0)  -  0,  -b  <  x  <  “  . 

Then  there  exists  a  constant  depending  only  upon  the  data  such  that 


<Pz> 


z(x,t)  <  /  rtl  |P  (x,tl-u(T) ,x) |dT 

0  X 


Therefore 


o(t)  a(t)  «  t  l5l 

/  v.(x,t)dx  <  /  z(x,t)dx  4  /  z(x,t)dx  <  K  /  -  dt  . 

0  -5(t)  -(t)  0 


The  estimate  of  I2  is  done  by  exploiting  again  the  methods  of  [2)  page  87  and 
dominating  the  integrand  by  z(x,t),  the  solution  of 

Lz  -  0,  (u(t)  <  x  <  "}  x  (0  <  t  <  T)  , 

z(x,0)  “0,  b  <  x  <  •  / 

z(a(t),t)  -  C^tt),  0  <  t  <  T  . 


It  gives 


22- 


t  Ijl 

J2  «  k14  /  7=  dT  • 

0  /t  -  T 


Proposition  5.2;  There  exist  a  constant  b  depending  upon  H,  b,  G^,  Gj,  T  such 
that 


Sg(t) 


l/  ufl(x,t)dx  -  /  ufl(x,t) dx |  <  B8  . 


Proof  of  Proposition  5 .2  s  For  t  >  0  fixed  we  have  Sg(t)  >  x^  .  then 

sQ(t)  x  x 

a  _  p  p 

|/  U0(x,t)dx  -  /  u0(x,t)dx|  <  /  |ug(x,t)  -  Ug(x,t )  | dx  + 

0  0  0  ” 

“a<t5  _ 

+  /  Ug(x,t)dx  “  J1  +  J2  • 


Standard  harrier  arguments  [7,8,9]  give 


sup  |Ug(x,t)  |  <  C1<b  +  C^T)  . 
(x,t)e0- 


On  the  other  hand  by  the  construction  of  b»(*),  the  distance  s-(t)  -  x  is  at  most 

O  w  p 

0^0,  Hence 

J2  4  k150  * 

To  estimate  the  Integrand  in  we  proceed  as  in  Lemma  4.1.  In  fact,  in  the  present 

situation  the  estimates  are  simpler  since  Ug  -  Ug  has  zero  flux  at  the  fixed  face 

x  •  0.  we  obtain  the  estimate  analogous  to  (4.1) 

|uB(x,t)  -  u„(x,t)|  <  K  sup  1 8 B ( T )  -  Sa(T)|  . 

®  9  16  0<T<t  9  9 


As  remarked  above  for  all  t  e  [0,T]  we  have 

S  g  ( t )  -  Sg(t)  <  , 

and  hence  the  proposition  Is  proved. 
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We  now  conclude  the  proof  of  lemma  5.1.  Notice  that  by  virtue  of  the  Lipschitz 


continuity  of  g(*,t)  and  (4.1)  we  have 


t  t,  t  |u(0,t)  -  Ufl<0,T)| 

/  |g(U(0,T)  ,T)  -  g(Ufl(0,t)  ,T)  |dT  <  G.T  '2  !  - - -  dT 


t  -  t 


1y  ,4,T 

<  G.T  /2c  /  .  dT  . 

0  /t  -  T 


Putting  together  the  various  estimates  so  obtained  we  see  that  for  all  o  <  t  <  T 
have 


t  iSl  t  161 

|s(t)  -  5g(t)|  <  K  /e  +  K  /  - -  dt  +  K  /  -  . T  dT  . 

0  /t  -  T  0  /t  -  T 


NOW 

is(t)  -  ;8(t)i  >  i s < t)  -  se<t)i  -  ise(t>  -  ie(t)i 

and 

l?l  -  sup  I s (p)  -  s-(p)|  <  ISI  +  eup  | s p ( P )  -  sfl(P)|  <  161  +  c,e  . 
0<P<T  0  1  0<P<T  0  0  T  1 


Consequently  the  above  implies 


151  <  K  /9  K  / 

t  20  21  i 


161 


dT 


The  proof  of  the  lemma  is  concluded  with  an  application  of  Gronwall's  inequality. 


we 
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6.  Modifications  of  the  scheme; 

Consider  the  sequence  of  approximating  problem  ( Pj )  introduced  in  Section  2.  The 
theorem  remains  true  if  we  modify  the  flux  condition  on  the  fixed  face  x  -  0  in  (P^) 
as  follows 

(6.i)  Retarded  flux. 

We  set  uj^(0,t)  »  g(h(0),t),  0  <  t  <  8  and  for  j  >  1  we  replace  the  flux 
condition  in  (Pj)  with 

ui(0,t)  -  g<uj-1(0,t  -  8),t),  (j  -  1)8  <  t  <  j8  . 

(6.ii)  Piecewise  constant  flux. 

The  modification  consists  in  freezing  the  flux  on  x  «  0  in  the  jth  rectangle 
at  its  value  at  the  lower  left  corner . 

Set  u^(0,t)  ■  g(h ( 0 ) ,0 ) ,  0  <  t  <  8  and  then  for  j  >  1 

u£(0,t)  -  glu^’tODj  -  1)0),  (j  -  1)8),  (j  -  1)0  <  t  <  j8  . 

The  proof  of  the  convergence  and  the  error  estimate  for  the  scheme  obtained  with  the 
above  modifications,  is  carried  out  in  essentially  the  same  way  as  indicated  in  the 
previous  sections.  Some  minor  modifications  are  needed  which  we  leave  to  the  reader, 
(cf.  [ID. 
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